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We show how the two-layer moist-convective rotating shallow water model (mcRSW), 
which proved to be a simple and robust tool for studying effects of moist convection 
on large-scale atmospheric motions, can be improved by including, in addition to the 
water vapour, precipitable water, and the effects of vaporisation, entrainment, and 
precipitation. Thus improved mcRSW becomes cloud-resolving. It is applied, as an 
illustration, to model the development of instabilities of tropical cyclone-like vortices. 
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l. Introduction 

Massive efforts have been undertaken in recent years in order to improve the quality of weather and climate modelling, and significant 
progress was achieved. Nevertheless, water vapour condensation and precipitations remain a weak point of weather forecasts, especially 
long-term ones. Thus, predictions of climate models significantly diverge in what concerns humidity and precipitations (Stevens and 
Bony 2013) . The complexity of thermodynamics of the moist air, which includes phase transitions and microphysics, is prohibitive. 
That is why the related processes are usually represented through simplified parameterisations in the general circulation models. 
However, the essentially non-linear, switch character of phase transitions poses specific problems in modelling the water cycle. 
Parametrisations of numerous physical processes in general circulation models often obscure the role of the water vapour cycle 
upon the large-scale atmospheric dynamics. The moist-convective rotating shallow water (mcRSW) model was proposed recently, 
precisely, in order to understand this role in rough but robust terms. The model is based on vertically averaged primitive equations 
with pseudo-height as vertical coordinate. Instead of proceeding by a direct averaging of the complete system of equations with full 
thermodynamics and microphysics, which necessitates a series of specific ad hoc hypotheses, a hybrid approach is used, consisting in 
combination of vertical averaging between pairs of isobaric surfaces and Lagrangian conservation of the moist enthalpy (Bouchut et al. 
2009; Lambaerts et al. 2011). Technically, convective fluxes, i.e. an extra vertical velocity across the material surfaces delimiting the 
shallow-water layers, are added to the standard RSW model, and are linked to condensation. For the latter a relaxation parametrisation 
in terms of the bulk moisture of the layer, of the type applied in general circulation models, is used. Thus obtained mcRSW model 
combines simplicity and fidelity of reproduction of the moist phenomena at large scales, and allows to use efficient numerical tools 
available for rotating shallow water equations. They also proved to be useful in understanding moist instabilities of atmospheric jets 
and vortices (Lambaerts et al. 2012; Lahaye and Zeitlin 2016; Rostami and Zeitlin 2017; Rostami et al. 2017). 

The mcRSW model, however, gives only the crudest representation of the moist convection. The water vapour can condense, but 
after that the liquid water is dropped off, so there are no co-existing phases and no inverse vaporisation phase transition in the model. 
Yet, it is rather simple to introduce precipitable water in the model, and link it to the water vapour through bulk condensation and 
vaporisation. At the same time, the convective fluxes present in mcRSW can be associated with entrainment of precipitable water, and 
its exchanges between the layers, adding more realism in representing the moist convection. Below, we will make these additions to the 
mcRSW model, and thus obtain an “improved” mcRSW, which we call imcRSW. We will illustrate the capabilities of the new model 
on the example of moist instabilities of hurricane-like vortices. Multi-layer modelling of tropical cyclones goes back to the pioneering 
paper Ooyama (1969), which had, however, a limited range due to the constraint of axisymmetry. Strictly barotropic models were 
also used, e.g. Guinn and Schubert (1993), as well as shallow water models with ad hoc parametrisations of latent heat release, e.g. 
Hendricks et al. (2014). The imcRSW model is a logical development of such approach. 
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Z 2 




Figure 1. Notations for the simplified two-layer scheme with mass flux across material surfaces. From Lambaerts et al. (2011), with permission of AIP. 



2. Derivation of the improved mcRSW 

2.1. Reminder on mcRSW and its derivation 


Let us recall the main ideas and the key points of derivation of the 2-layer mcRSW model. The starting point is the system of “dry” 
primitive equations with pseudo-height as vertical coordinate (Hoskins and Bretherton 1972). We recall that pseudo-height is the 
geopotential height for an atmosphere with an adiabatic lapse rate: z = zp ^1 — {p/po) R ^ Cp 'j, where zp = c p 9p/g, and the subscript 
0 indicates reference sea-level) values. Horizontal momentum and continuity equations are vertically averaged between two pairs of 
material surfaces zo, z\, and z \, Z 2 , where zp is at the ground, and z 2 is at the top. The pseudo-height z being directly related to 
pressure, the lower boundary is a “free surface” and the upper boundary is considered to be at a fixed pressure (“rigid lid”). The mean- 
field approximation is then applied, consisting, technically, in replacing averages of the products of dynamical variables by products of 
averages, which expresses the hypothesis of columnar motion. In the derivation of the “ordinary” RSW the fact that material surfaces 
z % , i = 0,1,2 are moving, by definition, with corresponding local vertical velocities allows to eliminate these latter. The main 
assumption of the mcRSW model is that there exist additional convective fluxes across Zj, such that 


wo 


dz 0 
dt ’ 


dz\ 

wi = — + Wi, 
dt 


W2 = + VV2, 

dt 


( 1 ) 


where Wi ,2 are contributions from the extra fluxes, whatever their origin, cf. Figure 1. The resulting continuity equations for the 
thicknesses of the layers h 2 = 22 — 21 , h\ = z\ — zq are modified in a physically transparent way, acquiring additional source and 
sink terms: 

d t hi + V • (hivi) = — Wi, 
dth,2 + V • (/ 12 D 2 ) = +Wi — VV 2 ■ 

The modified momentum equations contain the terms of the form WiV at the boundaries z,; of the layers. An additional assumption 
is, hence, necessary, in order to fix the value of the horizontal velocity at the interface. In the layered models the overall horizontal 
velocity, by construction, has the form v(z) = ViH(zi — z)H(z — where 7/(z) is Heaviside (step-) function. Assigning a 

value to velocity at z. t means assigning a value to the Heaviside function at zero, where it is not defined. This a well-known modelling 
problem, and any value between zero and one can be chosen, depending on the physics of the underlying system. In the present 
case this choice would reflect the processes in an intermediate buffer layer interpolating between the main layers, and replacing the 
sharp interface, if a vertically refined model is used. The “asymmetric” (non-centred) assignment 7/(0) = 1 was adopted in previous 
works. The “symmetric” (centred) assignment 7/(0) = 1 /2 will be adopted below. This choice does not affect qualitatively the previous 
results obtained with mcRSW, however it does affect the forcing terms in conservation laws. It corresponds to a choice of efficiency of 
momentum transport between the layers. In this way, the vertically averaged momentum equations become: 


dtv! + («i ■ V)«i + fk x V1 = —V0(zi) + g ^Vzi + 

d t v 2 + (v 2 ■ V)i> 2 + fk x v 2 = -V0(z 2 ) + Vz 2 + Wi + w 2 , 


Note that, whatever the assignment for Heaviside function, the total momentum of the two-layer system {z\ — zo)v± + (z 2 — z\)v 2 is 
locally conserved (modulo the Coriolis force terms). In what follows, will be assuming that W 2 = 0. 

The system is closed with the help of hydrostatic relations between geopotential and potential temperature, which are used to express 
the geopotential at the upper levels in terms of the lower-level one: 


<t>{zo) - zo) if20 < z < zi, 

0(Ao) TS'IjOi - zq) + g^(z - zi) ifzi < z < z 2 , 


The vertically integrated (bulk) humidity in each layer Qi = fqdz, i = 1,2, where q{x,y,z,t ) is specific humidity, measures 
the total water vapour content of the air column, which is locally conserved in the absence of phase transitions. Condensation introduces 
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a humidity sink: 


dtQi + V • (QiVi) = -Ci , i = 1,2. 


(5) 


In the regions of condensation (Ci > 0) specific moisture is saturated q(zi) = q s (zi) and the potential temperature 0(zi) + (L / c p )q s (zi) 
of an elementary air mass Wi dt dx dy , which is rising due to the latent heat release, is equal to the potential temperature of the upper 
layer Q i+1 \ 


Qi +1 


0(zi) + 



+ 


°p 


( 6 ) 


If the background stratification, at constant 9{zi) and constant q{zi), is stable 0i + 1 > by integrating the three-dimensional equation 
of moist-adiabatic processes 


d 

dt 



= 0. 


(7) 


we get 


Wi = ftCi, Pi 



1 

q(zi) 


> 0. 


( 8 ) 


In this way the extra vertical fluxes in (3), (2) are linked to condensation. For the system to be closed, condensation should be connected 
to moisture. This is done via the relaxation parametrisation, where the moisture relaxes with a characteristic time r c towards the 
saturation value Q s , if this threshold is crossed: 


Ci = Qi Q * H(Qi - Qi). (9) 

7"c 

Essentially nonlinear, switch character of the condensation process is reflected in this parameterisation, which poses no problem in 
finite-volume numerical scheme we are using below. For alternative, e.g. finite-difference schemes smoothing of the Heviside function 
could be used. In what follows we consider the two-layer model assuming that the upper layer is dry, and even with entrainment of 
water from the lower moist layer, water vapour in this layer is far from saturation, so the convective flux W 2 is negligible. In this way 
we get the mcRSW equations for such configuration: 

' dtV! + (wi • V)«! + fk X V! = - 5 V(hl + h 2 ) + 
d t v 2 + (v 2 ■ V)v 2 + fk x v 2 = -gS7{hi + sh 2 ) + Vl 2 ~” 2 fiC, 

< dth! + V ■ (hiwi) « -PC, , (10) 

dth 2 + V ■ (h 2 v 2 ) = +PC, 

. d t Q + V ■ (Qv 1 ) = -C, C= ^SlH(Q - Q s ) 

where s = 9 2 /9\ > 1 is the stratification parameter, vi = (m, i>i) and v 2 = ( u 2 , v 2 ) are the horizontal velocity fields in the lower and 
upper layer (counted from the bottom), with Uj zonal and ty meridional components, and hi , h 2 are the thicknesses of the layers, and 
we will be considering the Coriolis parameter / to be constant. 

As in the previous studies with mcRSW, we will not develop sophisticated parameterisations of the boundary layer and of fluxes 
across the lower boundary of the model. Such parameterisations exist in the literature (Schecter and Dunkerton 2009), and may be 
borrowed, if necessary. We will limit ourselves by the simplest version of the exchanges with the boundary layer, with a a source of 
bulk moisture in the lower layer due to surface evaporation E. The moisture budget thus becomes: 


d t Q + V- {Qv 1 ) = E — C 


The simplest parametrisations being used in the literature are the relaxational one 

e= QciRh{q-q), 

TE 


( 11 ) 


( 12 ) 


and the one where surface evaporation is proportional to the wind, which is plausible for the atmosphere over the oceanic surface: 

E oc |v|; (13) 

The two can be combined, in order to prevent the evaporation due to the wind to continue beyond the saturation: 

Q-Q 


Es = — -^|u| H(Q-Q). 

te 


(14) 


The typical evaporation relaxation time te is about one day in the atmosphere, to be compared with t c , which is about an hour. Thus 
te T c . Q can be taken equal, or close to Q s , as we are doing, but not necessarily, as it represents complex processes in the boundary 
layer, and can be, in turn, parametrised. 


(c) 2017 Royal Meteorological Society 
Prepared using qjrms4.cls 











112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 


Improved moist-convective rotating shallow water model 


5 


2.2. Improving the mcRSW model 


An obvious shortcoming of the mcRSW model presented above is that, although it includes condensation and related convective fluxes, 
the condensed water vapour disappears from the model. In this sense, condensation is equivalent to precipitation in the model. Yet, as is 
well-known, condensed water remains in the atmosphere in the form of clouds, and precipitation is switched only when water droplets 
reach a critical size. It is easy to include precipitable water in the model in the form of another advected quantity with a source due 
to condensation, and a sink due to vaporisation, the latter process having been neglected in the simplest version of mcRSW. We thus 
introduce a bulk amount of precipitable water, W(x, y , t), in the air column of a given layer. It obeys the following equation in each 
layer: 

d t W + V.(Wv) = +C - V, (15) 

where V denotes vaporisation. Vaporisation can be parametrised similarly to condensation: 

P= Qs ~ Q H(Q S -Q). (16) 

Tv 

Opposite to the condensation, vaporisation engenders cooling, and hence a downward convective flux, which can be related to the 
background stratification along the same lines as upward flux due to condensation: 


Wv 


-0*V, /3 * 


L* 

c v (o 2 -e o’ 


(17) 


where L* is the latent heat absorption coefficient, C v is specific heat of vaporisation. /3* is an order of magnitude smaller than f3. There 
is still no precipitation sink in (15). Such sink can be introduced, again as a relaxation with a relaxation time t p , and conditioned by 
some critical bulk amount of precipitable water in the column: 

W - Wrr 

P= — - —H{W~Wcr). (18) 

Tp 

The extra fluxes (17) due to cooling give rise to extra terms in mass and momentum equations of the model in each layer. Another 
important phenomenon, which is absent in the simplest version of mcRSW is the entrainment of liquid water by updrafts. This process 
can be modelled in a simple way as a sink in the lower-layer precipitable water equation, which is proportional, with some coefficient 
7 , to the updraft flux, and hence, to condensation, and provides a corresponding source of precipitable water in the upper-layer. 

Including the above-described modifications in the mcRSW models, and neglecting for simplicity 1) condensation and precipitations 
in the upper layer, by supposing that it remains far from saturation, 2 ) vaporisation in the lower layer, which is supposed to be close to 
saturation, we get the following system of equations: 


+ fz x vi = -gV(fti + h 2 ) + ( ® C V )( Vl 2 V2 ), 

<i ^ + fzxv 2 = -gV(h 1 + sh 2 ) + c h P V )( Vl 2 ^ 2 ), 
d t h 1 + V.{h 1 v 1 ) = -/3C + p*V, 

( d t h. 2 + V.(h 2 v 2 ) = +/3C -/3*V, ( 19) 

d t Wi + V.(Wiwi) = +(1 - 7) C - P, 
d t W 2 + V.(W 2 v 2 ) = +7C - V, 

OtQ i + V .(Qi^i) = — O + E, 

. 9tQ 2 + V.(Q 2 v 2 ) = V, 

where dj.../df = dt ... + (uj-V)..., i = 1,2. Here C is condensation in the lower layer considered to be close to saturation, W, is 
the bulk amount of precipitable water and Qi bulk humidity in each layer, 7 is the entrainment coefficient, V is vaporisation in the 
upper layer, considered as mostly dry. C, V, and P obey (9), (16), (18), respectively. Note that if the above-formulated hypotheses of 
mostly dry upper layer, and almost saturated lower layer are relaxed (or get inconsistent during simulations), the missing condensation, 
precipitation, and vaporisation in the corresponding layers can be easily restituted according to the same rules. 


2.3. Conser\’ation laws in the improved mcRSW model 


As was already said, the total momentum of the system is locally conserved in the absence of the Coriolis force (/ —> 0), as can be seen 
by adding the equations for the momentum density in the layers: 


(dt + vi.V)(hivi) + hiviV.vi + fz x (h\vi) = 

-g^Y ~ 9hlVh2 ~ ~ 


( 20 a) 


(dt + v 2 .V)(h 2 v 2 ) + h 2 v 2 V.v 2 + fz x (h 2 v 2 ) = 

-gs^ t f-gh 2 Vh 1 +('^Y^)(pC-fl*V) (2 ° b) 

The last term in each equation corresponds to a Rayleigh drag produced by vertical momentum exchanges due to convective fluxes. 

The total mass (thickness) h = hi + h 2 is also conserved, while the mass in each layer h\ 2 is not. However, we can construct a 
moist enthalpy in the lower layer 

h x - /3Qi - P*W 2 , (21) 


(c) 2017 Royal Meteorological Society 
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which is locally conserved, if entrainment and surface evaporation are absent: 

dtmi + V.(mi'Ui) = 0, i = 1,2. (22) 

The inclusion of precipitable water in the upper layer in (21) is necessary to compensate the downward mass flux due to vaporisation. 

The dry energy of the system E = f dxdy(e i + e 2 ) is conserved in the absence of diabatic effects, where the energy densities of the 
layers are: 

2 i 2 

j Vi . h i 

ei=hi-£ +g^~, 

< 

e 2 = h 2 ^ + ghih 2 + sg^. 

In the presence of condensation and vaporisation, the energy budget changes and the total energy density e = ei + e 2 is not locally 
conserved, acquiring a sink/source term: 


d t em -V ■ f e - (pC - p*V)g( 1 - s)h 2 , 


(23) 


where f e is the standard energy density flux in the two-layer model. For the total energy E = f dxdy e of the closed system we thus get 

d t E = (pC - P*V)g(s - 1) J dxdy h 2 . (24) 

For stable stratifications s > 1, the r.h.s. of this equation represents an increase (decrease) of potential energy due to upward (downward) 
convective fluxes due to condensation heating (vaporisation cooling). Note that with “asymmetric” assignment of Heaviside function 
at zero, an extra term corresponding to kinetic energy loss due to Rayleigh drag would appear in the energy budget, cf Lambaerts et al. 
( 2011 ). 

Potential vorticity (PV) is an important characteristics of the flow. In the presence of diabatic effects it ceases to be a Lagrangian 
invariant, and evolves in each layer as follows: 


di f Ci + / 

d t K hi 
d 2 / C 2 + / 

at ^ h 2 


Ci + / -, (pc — p*v) Z r (vi-v^pc^pvx 

hi ’' hi hi [ V 2 ’ fti J\ ’ 


C 2 + / x C 8C-P*V ) £ r ( vi-v 2 pc-p*v \ 

h 2 h 2 h 2 '[ V 2 ' h 2 ) 


(25a) 

(25b) 


where Q = i-(V x w,) = d x Vi — d y Ui ( i = 1, 2) is relative vorticity, and qi = (C,% + f)/hi is potential vorticity in each layer. One can 
construct a moist counterpart of potential vorticity in the lower layer with the help of the moist enthalpy (21), cf. Lambaerts et al. 
( 2011 ): 




Cl + / 

mi 


(26) 


The moist PV is conserved in the lower layer, modulo the Rayleigh drag effects: 


dl ( Cl+/ 

dt v mi 


= +z • 


V x 



pc-p*v y 

ml ). ■ 


(27) 


Note that the “asymmetric” assignment of the value of the step-function, which was discussed above, renders the moist PV in the lower 
layer conserved. 


3. Illustration: application of improved mcRSW model to moist instabilities of hurricane-like 
vortices 

3.1. Motivations 

We will illustrate the capabilities of the improved mcRSW, the imcRSW, on the example of moist instabilities of hurricane-like vortices. 
The mcRSW model, in its simplest one-layer version, captures well the salient properties of moist instabilities of such vortices, and 
clearly displays an important role of moisture in their development (Lahaye and Zeitlin 2016). Below we extend the analysis of Lahaye 
and Zeitlin (2016) to baroclinic tropical cyclones (TC), and use the imcRSW to check the role of new phenomena included in the 
model. Some questions which remained unanswered will be addressed, as well as new ones, possible to answer with the improved 
version of the model. In particular, we will investigate the influence of the size of TC (the radius of maximum wind) upon the structure 
of the most unstable mode, the role of vertical shear, and the evolution of inner and outer cloud bands at nonlinear stage of instability. 

3.2. Fitting velocity and vorticity distribution of the hurricanes 

We begin with building velocity and vorticity profiles of a typical TC within the two-layer model. An analytic form of the velocity 
profile is convenient both for the linear stability analysis, and for initialisations of the numerical simulations, so we construct a simple 
analytic fit with a minimal number of parameters. It has the form which is consistent with Mallen et al. (2005) where flight-level 
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Table 1. Parameters of the background vortices. BCW(S): weak(strong) baroclinic, BTW: weak “barotropic”, without vertical shear, l : the most unstable 
azimuthal mode 


config. 

l 

ei 

e 2 

Ql 

Ot2 

Pi 

Pi 

m\ 

m 2 

mi 

A) 2 

BCS 

3 

0.41 

0.49 

4.5 

4.5 

0.180 

0.178 

48 

47.5 

0.01 

0.0101 

BTW 

4 

0.4 

0.40 

2.25 

2.25 

0.25 

0.25 

14 

14 

0.1 

0.1 

BCW 

4 

0.4 

0.36 

2.25 

2.25 

0.25 

0.237 

14 

12.6 

0.1 

0.115 


observations were collected from Atlantic and eastern Pacific storms during 1977 — 2001: 


[(r-r 0 )“ i e- mi(r - ro) ' 3i ] 

e 2 , r t" , r ^ 

max (r — ro)“’e _mi (’ r o) /3s 
, m 0 r, r < r 0 . 


(28) 


Here i= 1,2 indicate the lower and upper layer, respectively, r is the non-dimensional distance from the center, t, measures the 
intensity of the velocity field, tq sets the non-dimensional distance of maximum wind from the centre, and other parameters allow to fit 
the shape of the distribution. A cubic Hermite interpolation across r = rg is made to prevent discontinuity in vorticity. Here and below 
we use a simple scaling where the distances are measured in units of barotropic deformation radius R d = \JgH //, and velocities are 
measured in units of \JgH , where H is the total thickness of the atmospheric column at rest. (Hence, the parameter e acquires a meaning 
of Froude number). Under this scaling the Rossby number of the vortex is proportional to the inverse of the non-dimensional radius of 
maximum wind ( RMW ). A useful property of this parametrisation is a possibility to tune the ascending or descending trends of the 
wind near and far from the velocity peak. Velocity is normalised in a way that the maximum velocity is equal to e. We suppose that 
velocity profile (28) corresponds to a stationary solution of “dry” equations of the model. Such solutions obey the cyclo-geostrophic 
balance in each layer: 

{v + f ) Vl = 9 Qr (i/l + H2) ’ (29a) 

(y + f ) V2 = 9 lr {Hl + aH2) ’ (29b) 

so the related Hi(r ) are obtained by integrating these equations using (28). The radial distribution of the relative vorticity in the vortex 
is given by (l/r)d [rV(r)] /dr. It should be emphasised that the radial gradient of the PV corresponding to the profile (28) has sign 
reversal, and hence the instability of the vortex is expected. Typical velocity and vorticity fields of an intense (category 3) vortex are 
presented in Figure 2. 
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0 12 3 

Normalized Radius (RMW Units) 



Figure 2. Normalised radial structure of azimuthal tangential wind with a fixed slope close to the centre (left panel) and relative vorticity (right panel) in both layers 
corresponding to the BCS vortex in Table 1. 


In what follows, we will be studying instabilities of thus constructed vortices, and their nonlinear saturation. The strategy will be 
the same as in (Lahaye and Zeitlin 2016): namely, we identify the unstable modes of the vortex by performing detailed linear stability 
analysis of the “dry” adiabatic system, with switched-off condensation and vaporisation, and then use the unstable modes to initialise 
numerical simulations of nonlinear saturation of the instability, by superimposing them, with small amplitude, onto the background 
vortex. We will give below the results of numerical simulations of developing instabilities for three typical configurations which are 
presented in Table 1: weak barotropic (BTW) and baroclinic (BCW), and strong baroclinic (BCS) vortices. 

3.3. Results of the linear stability analysis: the most unstable mode and its dependence on the radius of maximal wind 

By applying the standard linearisation procedure and considering small perturbations to the axisymmetric background flow, we 
determine eigenfrequencies and eigenmodes of the “dry” system linearised about the background vortex. The technicalities of such 
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analysis are the same as in Lahaye and Zeitlin (2015) extended to the two-layer configuration as in Rostami and Zeitlin (2017), and we 
pass directly to the results, which we will not give in detail either, limiting ourselves by what is necessary for numerical simulations in 
the next section. 

The most unstable mode with azimuthal wavenumber l = 3 of the BCS vortex is presented in Fig. 3. The unstable mode of Figure 3 




Figure 3. Upper row: Pressure and velocity fields in the x — y plane in the lower (left panel) and upper (right panel) layers corresponding to the most unstable mode with 
azimuthal wavenumber Z = 3 of the BCS vortex of Table 1. Lower row: left panel- corresponding divergence field of the most unstable mode, right panel - radial structure 
of three components of the most unstable mode: pressure anomaly 77, and radial (t>) and tangential (ii) components of velocity ; dashed (solid) lines: imaginary (real) part, 
thick (thin) lines correspond to upper (lower) layer. Note that the domain in the lower left panel is ss ten times larger than that of the upper panels. 


is clearly of mixed Rossby - inertia gravity wave type. Such unstable modes of hurricane-like vortices are well documented in literature, 
both in shallow-water models (Zhong and Zhang 2014) and full 3D models (Menelau et al. 2016). It should be stressed that at Rossby 
numbers which are about 40 and small Burger number, as can be inferred from the right panel of Fig. 2, the vortex Rossby wave part of 
the unstable mode, which is known since the work (Montgomery and Kallenbach 1997) and is clearly seen in the upper panels of Fig. 
3, is inevitably coupled to inertia-gravity wave field through Lighthill radiation mechanism, cf. Zeitlin (2008). 

The most unstable modes of BCW and BTW vortices have wavenumber l = 4. With our scaling the strength of the vortex is inversely 
proportional to its non-dimensional RMW, and thus the structure of the most unstable mode depends on RMW. Yet, as follows from 
Figure 4, the mode l = 4 is dominant through the wide range of RMW. In general, higher values of RMW correspond to higher 
azimuthal wavenumbers and lower growth rates. 

3.4. Nonlinear evolution of the instability 

We now use the unstable modes identified by the linear stability analysis to initialise numerical simulation of nonlinear evolution of the 
instability. We superimpose the unstable modes with weak amplitude (several per cent with respect to the background values) onto the 
vortex and trace the evolution of the system, as follows from numerical simulations with finite-volume well-balanced scheme developed 
for moist-convective RSW model (Bouchut et al. 2009). Numerical simulations with each of the vortex configurations of Table 1 were 
performed both in “dry” (M), with diabatic effects switched off, and moist-convective (MCEV) environments. 

The values of parameters controlling condensation, evaporation, vaporisation, and precipitation in MCEV environment are given in 
Table 2. It must be stressed that amount of precipitable water in each layer is highly sensitive to the values of parameters, especially to 
the intensity of surface evaporation. Condensation and precipitation time scales are chosen to be short, just few time steps At, while 
vaporisation and surface evaporation time-scales are much larger which is consistent with physical reality. Changing these parameters 
within the same order of magnitude does not lead to qualitative changes in the results. W cr is an adjustable parameter that controls 
precipitation and 7 controls entrainment of condensed water. The MCEV simulations were initialised with spatially uniform moisture 
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RMW (Rd Units) 

Figure 4. Dependence of the linear growth rates (in units of / _1 ) of the unstable modes with different azimuthal wavenumbers on the radius of maximum wind {RMW). 


Table 2. The values of parameters controlling condensation: t c , Q s , stratification: s = di/Bi, evaporation: T/,;, vaporisation: t v . and precipitation: W cr ,Tp 
At is the time-step of the code. 


T C 

Q s 

Q 

s 

te 

Tv 

Wcr 

Tp 

7 

5A t 

0.75 

~Q S 

1.5 

200Af «1 day 

10r c 

0.01 

3Af 

0.3 


content. We present below some outputs of the simulations, illustrating different aspects of moist vs “dry” evolution, and difference in 
behaviour of baroclinic and barotropic vortices. 

3.4.1. Evolution of potential vorticity 

We start by evolution of the PV field of the weak cyclone, as it is understandably slower than the evolution of the strong one, and 
different stages can be clearly distinguished. The evolution of potential vorticity in both layers during nonlinear saturation of the 
instability of the BCW vortex in MCEV environment is presented in Fig. 5, 6. The simulations show formation of a transient polygonal 
pattern inside the RMW at initial stages, with a symmetry of the most unstable linear mode. The patterns of this kind are frequently 
observed (Muramatsu 1986; Lewis and Hawkins 1982; Muramatsu 1986; Kuo et al. 1999). The polygon is further transformed into an 
annulus of high PV. Such annuli of elevated vorticity ( the so-called hollow PV towers (Hendricks and Schubert 2010) ) are found in both 
moist-convective and dry cases. It is worth mentioning that the growth of the primary unstable mode is accompanied by enhancement 
of outer gravity-wave field, as follows from the divergence field presented in Fig. 7. 

As follows from Figs. 5, 6 the polygon loses its shape at t « 17. At this time the modes with azimuthal wavenumbers l = 1,2 
are being produced by nonlinear interactions, and start to grow and interact with the polygonal eye-wall, which leads to symmetry 
loss by the core. A secondary, dipolar instability of the core thus develops, and gives rise to formation of an elliptical central vortex, 
corresponding to azimuthal mode l = 1, and of a pair of satellite vortices indicating the presence of l = 2 mode. The interaction of 
initial l = 4 mode with emerging l = 1 and l = 2 modes is accompanied by inertia-gravity wave (IGW) emission, and enhancement 
of water vapour condensation that will be discussed below. It should be emphasised that interaction between l = 2 mode and elliptical 
eye, of the kind we observe in simulations, was described in TC literature, e.g. Kuo et al. (1999), where reflectivity data from a Doppler 
radar were used to hypothesise that it was due to azimuthal propagation of l = 2 vortex Rossby waves around the eye-wall. 

Further nonlinear evolution consists in breakdown of the central ellipse with subsequent axisymmetrisation of the PV field, and its 
intensification at the center. This process characterises the evolution of both BTW (not shown) and BCW vortices, but is more efficient 
in the baroclinic case, as follows from Fig. 8. As seen in the Figure, the azimuthal velocity in the core region with r < 0.5RMW 
is subject to strong intensification. The exchanges of PV between the eye-wall and the eye, and intensification are known from the 
barotropic simulations (Montgomery and Kallenbach 1997; Schubert et al. 1999; Lahaye and Zeitlin 2016). As we see in Fig. 8, the 
intensification is enhanced by baroclinicity of the background vortex. This is confirmed by Fig. 9, and by Fig. 10, which illustrate the 
enhancement effect of both moist convection and baroclinicity upon palinstrophy, which is defined as 

V(t) = J f l VC.VCd.Td y, (30) 

in each layer, and which diagnoses the overall intensity of vorticity gradients. 

It is worth emphasising that, because of higher vorticity and smaller RMW, the axisymmetric steady state is achieved in the lower 
layer more rapidly than in the upper one in the case of baroclinic vortices. 

In the case of intense vortex, nonlinear evolution of the instability follows similar scenario, but is considerably accelerated, as follows 
from Fig. 11. 
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Figure 5. Nonlinear evolution of the most unstable l — 4 mode superimposed onto the background BCW vortex in MCEV environment, as seen in potential voiticity field 
in the lower layer. Formation of meso-vortices (zones of enhanced PV in the vorticity ring) is clearly seen, giving way to axisymmetrisation and monotonisation of the PV 
profile. Time is measured in units of / —1 . 


3.4.2. Spiral cloud bands 

Tropical cyclones exhibit specific cloud patterns. The new version of the model gives a possibility to follow clouds and precipitation, 
and it is interesting to test its capability to produce realistic cloud patterns. There are two types of cloud and rain bands associated with 
tropical cyclones, as reported in literature: the inner bands, that are situated close to the vortex core, within « 2 RMW, and the outer 
spiral bands located farther from the centre and having larger horizontal scales (Guinn and Schubert 1993; Wang 2009). Fig. 12 shows 
formation of inner and outer cloud bands, the latter having characteristic spiral form, during nonlinear evolution of the instability. Spiral 
cloud bands are related to inertia-gravity “tail” of the developing unstable mode. The link of spiral bands to inertia-gravity waves in 
“dry” RSW model of hurricanes was discussed in literature (Zhong et al. 2009). Here we see it in “cloud-resolving” imcRSW. It is to 
be stressed that amount of clouds strongly depends on the initial water vapour content. If it is closer to the saturation value, the amount 
of clouds and precipitation obviously increases and eventually covers the whole vortex. 

4. Conclusions and discussion 

Thus, we have shown that the moist-convective rotating shallow water model “augmented" by adding precipitable water, and 
relaxational parametrisations of related processes of vaporisation, precipitation, together with entrainment, is capable to capture some 
salient features of the evolution of instabilities of hurricane-like vortices in moist-convective environment, and allows to analyse the 
importance of moist processes on the life-cycle of these instabilities. There exist extended literature on the dynamics of the hurricanes 
eyewall, with tentative explanations in terms of transient internal gravity waves, which form spiral bands, cf. Lewis and Hawkins 
(1982) Willoughby (1978), Kurihara (1976), or alternative explanations Guinn and Schubert (1993) in terms of PV dynamics and 
vortex Rossby waves. Thus Schubert et al. (1999) obtained formation of polygonal eyewalls as a result of barotropic instability near the 
radius of maximum wind in a purely barotropic model, without gravity waves. A detailed analysis of instabilities of tropical cyclones 
was undertaken with a cloud-resolving model in Naylor and Shecter (2014), and showed that the results of Schubert et al. (1999) gave 
a useful first approximation for the eyewall instabilities. As was already mentioned in section 3.3, at high Rossby numbers the vortex 
Rossby wave motions are inevitably coupled to inertia-gravity waves, and our linear stability analysis confirms this fact. The mixed 
character of the wave perturbations of axisymmetric hurricane-like vortices was abundantly discussed in literature, e.g. Zhong et al. 
(2009). A detailed analysis of instabilities of hurricane-like vortices in continuously stratified fluid was given recently by Menelau et al. 
(2016), where it was shown that the inertia-gravity part of the unstable modes intensifies with increasing Froude number. The vortex 
profiles used above in section 3.2 have moderate Froude numbers, and the corresponding unstable modes have weak inertia-gravity 
tails. They are, however, sufficient to generate spiral cloud patterns, as we showed. The development of the instability of the eyewall 
proper at early stages (up to « 40/ _1 ) is only weakly influenced by moist convection, in accordance with findings of Naylor and 
Shecter (2014). This can be seen from comparison of the right and left panels of Fig. 9 and from Fig. 10. 

An advantage of our model, as compared to simple barotropic models, is its ability to capture both vorticity and inertia gravity waves 
dynamics. Another advantage is that the model includes the moist convection and couples it to dynamics and water vapour transport 
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Figure 6. Same as in Fig. 5, but for the upper layer. 
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Figure 7. Divergence field at t = 5[/ -1 ] in the lower (left panel) and upper (right panel) layers during the evolution of the instability of the BCW vortex in moist- 
convective and evaporating environment. Dashed line in the right panel represents the radius of maximum wind. Spiral structure, corresponding to inertia gravity wave 
“tail” of the unstable mode is seen to persist. 


278 and condensation in a self-consistent way, as compared e.g. with Hendricks et al. (2014), where the latent heat release due to the moist 

279 convection was introduced via ad hoc mass sink in rotating shallow water model. 
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Figure 8. Radial distribution of azimuthally averaged PV (left panel) and velocity (right panel) illustrating the role of the vertical shear in intensification of the vortex. 
Solid (dashed) lines represent the results of simulations in MCEV environment of the saturation of the instability of BCW (BTW) vortex in the lower layer. 




Figure 9. Effect of vertical shear on intensification of vorticity in the vortex core in environments without (M) and with {MCEV) moist convection and surface 
evaporation. The vorticity is normalised by its initial value. 




Figure 10. Evolution of palinstrophy in simulations with and without moist convection and surface evaporation of developing instability of the barotropic (left panel) and 
baroclinic (right panel) cyclones. 
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Figure 11. PV evolution during the saturation of / = 3 instability of the BCS vortex. Upper (lower) panel: lower (upper) layer. 
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Figure 12. Snapshots of precipitable water during the evolution of the BCS vortex displaying formation of inner and outer spiral bands in the lower layer. Thick dashed 
line represents RMW. Uniform distribution of initial Q is Q s — 0.005. 


280 Although we limited ourselves above by an application to tropical cyclones, the model can be used for analysis of various phenomena 

28 1 in mid-latitudes and tropics. The passage to the equatorial beta-plane is straightforward in the model, and it can be easily extended to the 

282 whole sphere. An important advantage of the model is that it allows for self-consistent inclusion of topography in the numerical scheme, 

283 giving a possibility to study a combination of moist and orographic effects. As was already mentioned, more realistic parametrisations 

284 of the boundary layer are available, and generalisations to three-layer versions are straightforward. 
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